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COMPLETELY MONOTONE FUNCTIONS AND SOME CLASSES 
OF FRACTIONAL EVOLUTION EQUATIONS 


EMILIA BAZHLEKOVA 


Abstract. The abstract Cauchy problem for the distributed order fractional 
evolution equation in the Caputo and in the Riemann-Liouville sense is studied 
for operators generating a strongly continuous one-parameter semigroup on a 
Banach space. Continuous as well as discrete distribution of fractional time- 
derivatives of order less than one are considered. The problem is reformulated 
as an abstract Volt err a integral equation. It is proven that its kernel satisfy 
certain complete monotonicity properties. Based on these properties, the well- 
posedness of the problem is established and a series expansion of the solution is 
obtained. In case of ordered Banach space this representation implies positivity 
of the solution operator. In addition, a subordination formula is obtained. 

Keywords: completely monotone function, Bernstein function, fractional evo¬ 
lution equation, distributed order fractional derivative, Mittag-Leffler function 


1. Introduction 

The evolution equations of fractional order are extensively used for modeling of 
materials and processes with memory. Recently, in the attempt to find more ade¬ 
quate models, equations involving discrete or continuous distribution of fractional 
derivatives are introduced. 

In this paper, we consider the fractional evolution equation of distributed order 
in the following two alternative forms: 

l 

/.i(/3) c ufw(t) d(3 = Au(t ), t > 0, (1.1) 

and 

u'(t)= f n(j3)Dt Au(t) d/3, t > 0, (1.2) 

Jo 

where C D^ and Df are the fractional time-derivatives in the Caputo and in the 
Riemann-Liouville sense, respectively, and A is a closed linear unbounded operator 
densely defined in a Banach space X. The initial condition it(0) = a £ X is 
prescribed. For the weight function /r two cases are considered: 

• discrete distribution 

m 

M/3) = 6(P - a) +'^2b j S(P - aj), (1.3) 

j'=i 

where 1 > a > or... > a m > 0, bj > 0, j = 1, m > 0, and S is the 

Dirac delta function; 
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• continuous distribution 

/*eC[0,i], n(J3) >o, [0,1], (1.4) 

and /i(/3) / 0 on a set of a positive measure. 

In the case of discrete distribution, equations ed and m are reduced to the 
multi-term time-fractional equations 

m 

c D^u(t) +’ s ^bj C Dt 3 u(t) = Au(t), t> 0, (1.5) 

3 =1 

and 

m 

u'(t)=D?Au(t) + Y l b 3 D?Au(t), t> 0, (1.6) 

i=i 

respectively. Note that if m = 0 (single-term fractional evolution equation) problem 
ED is equivalent to ed with a replaced by 1 — a. However, in general, similar 
equivalence does not hold for equations ED and ED- 

Problems ED and ED have found numerous applications in the modeling 
of anomalous relaxation and diffusion phenomena, see e.g. 02 m [29j. This 
leads to a growing interest to such problems in the last years. Concerning the 
equation in the Caputo sense ED, we refer to 0 E M E US ED for the study 
of its multi-term variant ED and to 0 E E HU E for the case of continuous 
distribution. Concerning the equation in the Riemann-Liouville sense ED , two- 
and three-term variants of problem ED are considered in E HI HU and the 
continuous distribution case of ED is discussed in 22s : 22U22]- Although problem 
ED is already studied for various classes of operators A: the Laplace operator in 
different settings EE, the Riesz space-fractional derivative E, second order 
symmetric uniformly elliptic operators EHHH, general operators in a Banach 
space setting E, etc., to the best of the author’s knowledge, the mathematical 
study of problem ED is still quite limited. 

In this paper, it is assumed that the operator A is a generator of a Co-semigroup 
(see e.g. 0),i.e. that the classical abstract Cauchy problem 

u'(t) = Au(t ), t > 0; u( 0) = a £ X, (1-7) 

is well-posed. Reformulating problems ED and ED as abstract Volterra integral 
equations, we propose a unified approach to their study. We prove that the scalar 
kernels of the corresponding integral equations have certain complete monotonicity 
properties and derive useful consequences for the original equations ED and ED 
based mainly on these properties. 

Recall that a function / : (0, oo) —> R is called completely monotone if it is of 
class C°° and 

(—1 ) n /(™)(t) > 0, for all i > 0, n = 0,1,... 

Completely monotone functions appear naturally in the models of relaxation and 
diffusion processes. Such are the exponential and the Mittag-Leffler functions of 
negative argument, which are obtained as solutions of the classical and the single¬ 
term fractional relaxation equations. For details and references on the complete 
monotonicity property see e.g. (26] and EI3, Chapter 4. Concerning Mittag-Leffler- 
type functions we refer to iii- 

The present paper is organized as follows. Section 2 contains preliminaries. In 
Section 3 problems ED and ED are rewritten in equivalent form as abstract 
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Volterra integral equations and the properties of their scalar kernels are studied. 
In Section 4 the wcll-posedness of the problems is proven and series expansion 
of the solutions is obtained in terms of the resolvent operator of A. In addition, 
the positivity of the solution operators is discussed. In Section 5, a subordination 
formula is obtained. 


2. Preliminaries 


The sets of positive integers, real and complex numbers are denoted as usual by 
N, R, and C, respectively, and No = NU{0}, R+ = [0, oo). Denote by Eg the sector 

Eg := {s e C; 0, | args| < d}. 


Denote by * the Laplace convolution: 

{h*h){t)=( fi(t — t )/ 2 ( t ) dr. 

Jo 

Let J t 7 be the fractional Riemann-Liouville integral 

J * f{t) = TW)J 0 ( i -^V(r)dr={L^}*/(t), 7 >0, 

where T(-) is the Gamma function. 

The Caputo and the Riemann-Liouville fractional derivatives of order P G [0,1], 
c Dt and Z)f, are defined by C D® = = /, C D\ = D] = d/dt , and 

c D^ = Jl~ p Dl = P G (0,1). 


Application of the Laplace transform 

pO O 

C{f(t)}(s) = f(s)= e~ st f(t) dt 
Jo 

to the operators of fractional integration and differentiation gives m 

7 > 0, 

C{ c D P t f}(s) = aPf(s) - s /3_1 /(0+), 0 G (0,1), 
£{D?f}{a) = a p f(a) - (J}-*f){0+), P € ( 0 , 1 ). 
Denote as usual by E a> p(-) the Mittag-Leffler function 


E a A z ) = 


fc=o 


T(ak + PY 


a, p,z G C, Jia > 0. 


Recall the Laplace transform pairs (IRa >0, /iGl, t > 0) 


4-Ct—l 


T(a) 


= s~ a , C{t p - 1 E atP (- f i1? l )} = 




( 2 . 1 ) 


( 2 . 2 ) 


The characterization of completely monotone functions is given by the Bern¬ 
stein’s theorem (see e.g. J5j) which states that a function / : (0,oo) —> R is 
completely monotone if and only if it can be represented as the Laplace transform 
of a nonnegative measure. 

A C°° function / : (0, oo) —>■ R is called a Bernstein function if it is nonnegative 
and its first derivative f'(t) is a completely monotone function. 
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The classes of completely monotone functions and Bernstein functions will be 
denoted by CMT and BT. Some properties of these classes of functions, which will 
be used further, are summarized in the next proposition. 

Proposition 2.1. The following properties are satisfied: 

(a) The class CMT is closed under pointwise addition and multiplication; 

(b) If f € CMT and <p £ BT, then the composite function f(ip) £ CMT; 

(c) If ip £ BT, then <p{s)/s £ CMT; 

(d) Let f £ L} oc (WL+) be a nonnegative and nonincreasing function, such that 
lim t _> +00 f(t) = 0. Then tp(s) = sf(s) £ BT; 

(e) If f G T; 1 oc (-^+) an d f ® s completely monotone, then f(s) admits analytic ex¬ 
tension to E,r and |arg/(s)| < |args|, s £ E T . 

Proof. Properties (a) and (b) follow by application of the product and the chain 
rule of differentiation (a detailed proof of (b) can be found in [26]). For (c) and (d) 
see the proof of Proposition 4.3 in EZj. For (e) see m, Example 2.2. □ 

Let X be a Banach space with norm ||.||. Let A be a closed linear unbounded 
operator in X with dense domain D(A), equipped with the graph norm ||. ||>i, 
||x||a := ||x|| + ||Ax||. Denote by g{A) the resolvent set of A and by R(s,A) the 
resolvent operator of A: R(s,A ) = (s — A ) _1 , s £ g(A). 

Next we recall some definitions and basic theorems, given in EH, concerning the 
abstract Volterra integral equation 

u(t) = a + f k(t — t)Au(t) dr, t > 0; a £ X, (2.3) 

Jo 

with a scalar kernel k £ Lj oc (R + ). 

A function u £ C(M + ; X) is called a strong solution of (12.311 if u £ C(R + ; D(A)) 
and (12.31) holds on R + . 

Problem (12.31) is said to be well-posed if for each a £ D{A) there is a unique 
strong solution u(t;a) of (12.31) and a n £ D(A), a n —> 0 imply u(t;a n ) -A 0 in X, 
uniformly on compact intervals. 

For a well-posed problem the solution operator S(t) is defined as usual by 
S(t)a = u(t ; a), a £ D(A), t > 0. 

In this paper it is assumed that A is a generator of a bounded Co-semigroup T(t), 
||T(£)|| < M, t > 0. The Hille-Yosida theorem (see e.g. [S]) then implies that 
(0, oo) C g(A) and 

|| R{s, A)" || < M/s n , s > 0, n £ N. (2.4) 

Suppose / 0 °° e~ st \k(t)\dt < oo for s > 0 and k(s) ^ 0, 1 /k(s) £ g(A) for s > 0. 
Then the Laplace transform of the solution operator S(t) of problem (12.31) 

pOO 

H(s) = / e~ st S(t)dt , 3fo>0 

Jo 

is given by 

H i s ) = ^y-R{g(s),A), g(s) = l/k(s). (2.5) 

The Generation Theorem ([27], Theorem 1.3) states that problem (12.31) is well-posed 
with solution operator S(t) satisfying ||S(t)|| < M, t > 0, if and only if 


for all s > 0, n £ N 0 . 


( 2 . 6 ) 
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3. Integral reformulation and properties of the kernels 

We reformulate problems EH and (11.211 as Volterra integral equations of the 
form (12.31) with appropriate kernels k(t). By applying (formally) the Laplace trans¬ 
form and, by the use of EH, it follows for the solution of (12.31) 

u(s) = -(1 — k(s)A)~ 1 a (3.1) 

s 

and for the solutions of problems EH and EH, respectively 

-Aj a, (3.2) 

where 

h(s) = ( n{(3)sP dp. (3.3) 

Jo 

Note that in the discrete distribution case (11.31) h(s) admits the representation 

m 

h(s) = s a + J2 b jS aj . (3.4) 

i=i 

Comparing (13.21) to (13.11) . it follows for the kernels fci(f) and A 12 (^), corresponding 
to problems EH and EH, respectively: 

ki(s) = (h(s )) -1 , k 2 {s) = h(s)/s. (3.5) 

Representation (12.51) suggests that it is convenient to define also the functions: 

g l {s) = l/ki{s), i = 1,2, 

that is 

5 i(s) = h(s), g 2 {s) = s/h(s), (3.6) 

where h(s) is defined in (13.31) . Some useful complete monotonicity properties of the 
functions ki(t) and gi(s), i = 1 , 2 , are proven in the next theorem. 

Theorem 3.1. Let p(p) be either of the form 11.3ft or of the form EH with the 
additional assumptions p G C 3 [0,1], p(l) ^ 0, and p(0) / 0 or p(p) = ap v as 
P —>■ 0, where a,v > 0. Then the functions ki{t) and gi{s), i = 1,2, have the 
following properties: 

(a) k, G and lim^ +00 kiit) = 0; 

(b) ki{t ) G CM.T for t > 0; 

(c) ki* k 2 = 1; 

(d) gi(s ) G BT for s > 0; 

(e) gi(s)/s G CM.T for s > 0; 

(f) 9i{ s ) admits analytic extension to and for any sGS, 

|arg 5 l (s)| < | args|. 

In the multi-term case (3.4\ ) a stronger inequality holds: 

I arg^(s)| < a\ args|, sG £*•. 

Let us first consider some particular cases. Applying EH , EH , EH and (liTfil) . 
it follows in the single-term case ( (13.41) with m = 0): 

J-QL— 1 4- — CK 

ki{t) = tw r, k 2 (t) = —--r, 51 (s) = s a , g 2 (s) = s 1_ “, 

1 (a) 1 (1 — a) 


u(s) = 1^1 (h(s) — A) 1 a, u{s) = — 
s s 


1 - 


h(s 
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and in the double-term case (' (13.41) with m = 1): 


h(t) = r- 1 £ a _ ai , a (-M“- ai ), k 2 (t) = 


gi(s) = s a + b 1 s ai , g 2 (s ) 


s 

s a + bis ai 


t -a ' ^ t - ai 

r(l -a) + 6 l r(l-ai)’ 


Thus, in the single-term case Theorem 3.1 is straightforward. In the double term 
case, statements (a) and (b) are trivial for k 2 ] for k\ they follow from the fact that 
the Mittag-Leffler function E a< p(—x) G CMP for x > 0, 0 < a < 1, /3 > a (see e.g. 
[25] ) and Proposition 2.1(a) and (b). On the other hand, properties (d) and (e) are 
trivial for g i, but need more elaborate proof for g 2 . 

Finally, consider the case of continuous distribution in its simplest form: constant 
weight function /r(/3) = 1. Then (13.31) implies (taking s 13 = e^ logs ) 


5i 00 


s — 1 
logs ’ 


52 (s) 


slogs 
s — 1 


Based on these explicit representations, only the positivity of the functions <7i(s) 
for s > 0 is evident, however properties (d) and (e) are not easily recognized. 

Now, we proceed with the proof of Theorem 3.1. 


Proof. Let us start with the kernel k 2 (t). Application of the inverse Laplace trans¬ 
form to k 2 (s) = h(s)/s , see (13.51) . implies by the use of (12.21) : 

k 2 (t) = m d@. (3.7) 

In the case of discrete distribution, (13.71) and (13.41) imply 


k 2 {t) = 


t~ a 

r(i — ol) 



r(i -<*,)' 


Therefore 

k 2 (t) ~ t ~ a , t —> 0; k 2 {t) ~ t~ am , t —> oo, 
and thus (a) is satisfied for this kernel. 

In the case of continuous distribution it is proven in nu, Proposition 2.1, that 


k 2 {t) 


1 

i(bifF 


t 0. 


Therefore, it is integrable at t = 0 (note that the singularity at t = 0 is quite 
strong). Moreover, since T(1 — /3) > 1 for /3 G [0,1], (13.71) implies for t > 1 

0 <k 2 {t)< sup |/i(/3)| [ —~~- 

/3e[o,i] Jo tlogt 


and thus k 2 {t) —> 0 as t —> oo. Complete monotonicity of k 2 {t) follows directly by 
differentiation of (13.71) . In this way, (a) and (b) are proven for the kernel k 2 (t) in 
both discrete and continuous case. 

Consider now the kernel hi. The identity fci(s) = l/h(s), see (13.51) . implies the 
following representation for this kernel as an inverse Laplace integral: 


ki(t) 


1 

27ri 


/»7+ioo 



ds, 


7 > 0 . 


(3.8) 
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Consider first the discrete case in which h(s) is defined by (13.411 . The function h(s) 
has no zeros in E T . Indeed, for s = re 1 ^ , with r > 0, </> € (— n,n), 

m m 

S{s“ + E bjS aj } = r“ sin cu/> + E sin aj4> 7 ^ 0 , 

7=1 7=1 

since sina</> and sin aj(f> have the same sign and bj > 0. Then the function under 
the integral sign in (13.81) is analytic in T, n and we can bend the integration contour 
to the contour T p ^ defined by 


T p-s = T P ,e u r °,e u r o.e> P > ?r/2 < 6 < tt, 


pfi 


(3.9) 


where 


T± e = {re ±ie : r > p) , T° e = {pj* : < 6} , 


and Tp.e is oriented in the direction of growth of arg s. Hence 




1 


e — 

c o 

r„.o s 


■EZibB 


■ ds. 


(3.10) 


The integral over T° g is a function from C°°[0, 00 ). Take p = R so large that 

m m 

k“ + E 6 / ] 1 > ^ M a A \s\>R. 

7=1 7=1 

Then, noting that cosd < 0 for 7 r /2 < 0 < n, it follows 

1 


Ir~ 1 ir+ 

1 R,( 




■ ds 


S 3 


/•OO 

<C e rtcose r - a dr < C ^-l_ ( 3-11 ) 

Jr 


Therefore, k\{t) ~ i“ 1 for t —► 0 and thus it has an integrable singularity at f = 0. 
Since in the discrete case k\(s) ~ s~ am as s —> 0, Karamata-Fcller Tauberian 
theorem ( 0 , Chapter XIII) implies k±(t) ~ t —> 00 . Thus (a) is proven 

for the discrete variant of ki(t). To prove its complete nronotonicity we take p —> 0 
and 9 —> n in (13.101) . Since 


r° s 

1 p,e 




■ ds 


< C 


e ptcos V“ m #, 


(3.12) 


the integral over T p g vanishes when p —7 0. Therefore, only the contributions of 
the integrals over T^ g remain in (13.101) . implying 


/•OO 

ki(t) = / e~ rt K(r)dr , 

Jo 


(3.13) 


where 


K(r) =-3 1 


Simplifying this expression, we get 


1 


K(r) = — 


Z?= iM c 

B{r) 


tt (H(r)) 2 + (B(r)Y 


where 


A{r) = r" cos an + E bj r0lj cos ajn, B(r) = r" sin an + E ^j r0tj si n a 7 7r ) 
7=1 7=1 
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and thus K(r) > 0 for r > 0. This together with representation (13.1311 implies that 
ki(t) G CMT. 

In the case of continuous distribution it is proven in m, Proposition 3.1, that 
for small values of t 

ki(t) < Clogi, 

therefore this kernel has integrable singularity at t — > 0. Further, by m , Proposi¬ 
tion 2.2, (ii) and (iii), 

/ i\ A +! 

fci(s) ~ ( lo S — ) > s->0, 


where 


A = 


0 if/z(0)^0, 
v if g((3) = aft” as fi 0. 


Applying again Karamata-Feller Tauberian theorem m, Chapter XIII) it follows 

, ^ (logt) A , , 
ki(t) ~ ---, t-S' oo, 


and thus k\{t) —> 0 as t —>• oo. Complete monotonicity of ki(t) in the case of 
continuous distribution is proven in Propositions 3.1. In this way the proof of 
properties (a) and (b) is completed for all cases. 

According to (13.511 

ki(s)k 2 (s) = 1/s 

and taking the inverse Laplace transform of this identity we derive (c). 

Further, gi(s) G BT and gi{s)/s G CMT by direct differentiation (see (j3.6|l and 
(13.311 b To prove that <72(s) G BT we use its representation (see (13.511 and (13.61) 1 

g 2 (s) = ski(s). 


In view of the properties of ki(t) proven in (a) and (b) we can apply Proposition 
2.1(d) which implies that g 2 (s) is a Bernstein function. Then by Proposition 2.1(c) 
g 2 (s)/s G CMT. In this way statements (d) and (e) are proven. 

The first part of (f) follows from the complete monotonicity of k\(t) and k 2 (t), 
applying Proposition 2.1(e). The stronger mapping property of gt in the multi- 
term case follows by applying the inequalities | arg(s /3 )| = /3| args|, | arg(si + S 2 )| < 
max{| argsi|, | args 2 |} and arg(s _1 ) =args. □ 


Let us now prove that problem HU is equivalent to the abstract integral equa¬ 
tion (12.311 with kernel k = k± and problem (11.211 is equivalent to (12.311 with kernel 
k = k 2 . Following ideas from nu, we can rewrite the Caputo differential operator 
of distributed order 

C B^= [\(P ) c D p t dp 
Jo 

in the form 

c O^/ = k 2 * f = (*2 * /)' - k 2 (t)f( 0), (3.14) 

for functions / for which the expressions are well defined. Consider also the integral 
operator of distributed order for / G L 1 

4 M) / = h * /. 


( 3 . 15 ) 
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The functions k\ and k 2 in (13.151) and (13.141) are exactly the kernels, corresponding 
to problems (11.11) and (11.21) . Then 

C D = jW c D^ ) / = /-/(0). (3.16) 

Indeed, (13.141) . (I3.15[) and Theorem 3.1.(c) imply 

= (k 2 * (h * /))' - fc 2 (*)(*! * /)(0) = ((fca * h) * /)' = (1 * /)' = /, 

jWc B W / = ki * (fc2 * = ( fcl * fc 2 ) * /' = 1 * /' = f _ /(0). 


Using identities (13.161) . we derive the Volterra integral equation (12.31) with kernel 
k = ki by applying operator to both sides of equation HU- Conversely, 
applying to both sides of (|2.3D with k = k\ , we get back equation (11.11) . The 

initial condition u(0) = a is recovered by taking t = 0 in (12.31) . 

Concerning the Riemann-Liouville differential operator of distributed order 

]©£ m) = t dp, 

Jo 

it can be represented in the form 

D ^ } / = (*2 * /)'- 


This implies J i I ID)j ,i ' 1 / = k 2 * /■ Then, the equivalence of problem m and the 
abstract integral equation (12.31) with kernel k = k 2 can be established in a similar 
way as above, but applying the classical integration and differentiation operators 
J} and Dj instead of and c O^. 


4. WELL-POSEDNESS, REPRESENTATION FORMULA AND POSITIVITY 

The following definition is based on the equivalence proven at the end of the 
previous section. 

Definition 4.1. Problem Hl.l\) . resp. mu, is said to be well-posed if the Volterra 
integral equation 12.3 1) with kernel k = k\, resp. k = k 2 , is well-posed. The solution 
operator of the corresponding integral equation \2.3\) is called a solution operator of 
problem \1.1\) . resp. El- 

Now we formulate the main result of this paper. 

Theorem 4.2. Suppose that the conditions of Theorem 3.1. on the weight function 
/u(/3) are satisfied. Let A be a generator of a bounded Co-semigroup T(t), such that 
||T(f)|| < M, t > 0. Then problems if 1.1]) and HI.21) are well-posed and their solution 
operators S(t) satisfy ||S'(t) || < M, t > 0. Moreover , the solutions of 11.1)1 and \1.2 1) 
have the representation 

1 n k 

u(t)= lim — (n/t) n+ls ^^b n , k}P (n/f)(R(g{n/t) ,A)) P+1 a, (4.1) 

n—too Ti\ z —' z — y 

k —0 V— 1 

where the convergence is uniform on bounded intervals of t > 0. The functions 
b n ,k,p(s) are nonnegative for s > 0 and are defined by 



b n ,k, P (s) = (-l) n+p 


a k , P (s)p\, s > 0, 


(4.2) 
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where the functions a k , P (s) are given by the recurrence relation 

a k +i, p {s) = a k , P -i(s)g'(s) + a kp (s), 1 < p < k + 1, k > 1, (4.3) 

a k ,o & k ,k +1 = 0? ^i,i(^) 9 (^)* 

Here g(s) = gi(s) in case of problem \1.1\ ) and g(s) = < 72 ( 5 ) in case of problem 

m- 


Proof. First we prove (EH), where H(s) is defined by (12.51) with g(s) = gi(s), i = 
1, 2. Let us express H( n \s) in terms of powers of 

w(s) = R(g(s), A). 


Note that by Theorem 3.1 if s > 0 then g(s) > 0 and thus g(s) G 9 (A), i.e. the 
resolvent operator R(g(s),A) is well defined. By the Leibniz rule it follows 


h^\ S ) = j2 


k—0 


g(s) 


( n—k) 


wW(s). 


The formula for the fc-th derivative of a composite function (see [30]) gives 
w [k \s) ='^2a k . p {s)(-l) p p\{R{g(s),A)) p+1 , 

P= 1 

where the functions dk, p (s) are defined by (POl) . 

We will prove inductively that for any k > 1 and 1 < p < k 

{-l) k+p a k , p (s ) G CMT. 


(4.4) 


(4.5) 


(4.6) 


For k = p = 1 this is fulfilled since aip(s) = g'(s) and g'(s) G CMJ- by Theorem 
3.1(d). Further, 02,1 = g ", 02,2 = ( g ’) 2 and the assertion (|4.6I) holds for these 
functions applying Theorem 3.1.(d) and Proposition 2.1.(a). Now fix some ko > 2 
and suppose that (14.61) holds for all k < ko, 1 < p < k. Then, (14.31) implies that (14.61) 
is satisfied for k = ko + 1, 1 < p < ko, since (—l) fc ° +p+1 afc 0lP -i(s)( 7 , (s) G CMP as 
a product of two completely monotone functions and (—l) fco+p+1 a' fco p (s) G CMP 
by (]T 6 J) . In addition, by (| 1 ■ • 1 P , a k ^~r. a,q —i— 1 — cA:,,. k,, g and it is completely monotone 
since a ko , ko G CMP and g’ G CMP. In this way the proof of (14.61) is completed. 

In particular, (14.61) implies 

(-l) fc+p a fc , p (s) >0, s > 0. (4.7) 

On the other hand, by Theorem 3.1(e) g{s)/s G CMP, i.e. 

(-I)""* ( 9 (s)/s) {n - k) > 0, s > 0. (4.8) 

Inserting (14.51) in (14.41) we get 

n k 

(-1 ) n H (n \s) ( 4 . 9 ) 

k—0 p—1 

where the functions b n ,k, P {s ) are defined in (|4.2D . Moreover, inserting m and 
EH in EH, it follows 

b n ,k,p(s ) >0, s > 0. (4-10) 

In addition, let us note that in the trivial case A = 0 (14.91) implies the identity 

n k 

(-ins- 1 )™ = E v*»G?(s)r Cp+1) • 

k—0 p—1 


(4.11) 
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Now, applying successively (14.101) . (12.41) and (14.111) we obtain from (14.91) 

n k 

\\H {n \s)\\ 

k —0 p= 1 
n k 

< mJ2J2<^ s ^9(s))- { p+1) 

k =0 p—1 

= M(-l) n (s- 1 )W=Mn!r (n+1) , s > 0. 

Therefore, conditions (12.6|) are satisfied and by the Generation Theorem problem 
m is well-posed with bounded solution operators S(t), satisfying ||5(t)|| < M, 
t > 0. Therefore, by Definition 4.1 this holds for problems CD and m- 

Further, we use the Post-Widder inversion formula for the Laplace transform 
(see e.g. my Let u(t), t > 0, be a X valued continuous function, such that 
u(t) = 0(e 7t ) as t —>• oo for some real 7. Then 


u(t ) = lim 

n—>• 00 


(-1)" 

n\ 


(?r (£) (?) 


uniformly on compact subsets of (0, 00). 

Since u{t) = S(t)a is a continuous and bounded function for t > 0, the Post- 
Widder inversion theorem can be applied and gives the representation CD- □ 


For an analogous result concerning the single-term fractional evolution equation 
see [2], Corollary 2.10. 

The representation CD is a generalization of the exponential representation for 
the solution of the classical Cauchy problem CZD 

u(t) = lim (i -^4] a. 

n— J-oo y n ) 

The positivity of the coefficients b n ^^ p in representation CD has a useful direct 
consequence: it implies the positivity of the solution operator. This holds, however, 
only in spaces in which the notion of positivity is well-defined, i.e. in ordered Banach 
spaces (for a simple introduction see e.g. jSJ, p. 353). 

Suppose X is an ordered Banach space. For example, such are the spaces of type 
L P (Q) or Co (ft) for some fi G R d , d G N, with the canonical ordering: a function 
a G A' is positive (in symbols: a > 0) if a(x ) > 0 for (almost) all ie!l. 

A solution operator S(t) in an ordered Banach space X is called positive if a > 0 
implies S(t)a > 0 for any t > 0. 

In other words, positivity of a solution operator means that positivity of the 
initial condition is preserved in time. Next we prove that this is satisfied for the 
considered problems CD and m if the operator A generates a positive Co¬ 
semigroup. 


Corollary 4.3. Let X be an ordered Banach space. Assume the conditions of 
Theorem f.2. are satisfied and the solution operator T(t ) of the classical Cauchy 
problem jl.7\ ) is positive. Then the solution operators S(t) of problems \1.1\) and 
m are positive. 


Proof. Since 


R{s,A) 


e~ st T{f) dt , 


s > 0, 
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the positivity of the Co-semigroup T(t) imply that the resolvent operator R(s, A) is 
positive: if a £ X and a > 0, then R(s,A)a > 0, s > 0. Therefore R(g(s),A)a > 0 
for all s > 0. This together with the positivity of the coefficients (14.21) in the 
representation formula m implies the positivity of S(t). □ 

Positivity related to problem ED with operator A = A on R d , d > 1, is estab¬ 
lished in jH [29l M ill) by proving the positivity of the fundamental solution. In 
the case of a second order elliptic operator A on a bounded domain, positivity is 
implied by a corresponding maximum principle, see e.g. |21| . 


5. Subordination formula 


In the previous section we proved that well-posedness of the classical Cauchy 
problem ED implies well-posedness of problems ED and ED- For completeness, 
here a subordination formula is given, relating the solution operator S(t) of the 
fractional evolution equation ED , resp. ED. and the solution operator T[t) of 
the classical Cauchy problem ED- 


Theorem 5.1. Assume the conditions of Theorem 3.1. on the weight function 
p,{fJ) are satisfied and let A he a generator of a bounded Co-semigroup T(t). Then 
the solution operator S ( t ) of problem ili.il) . resp. ED. satisfies the subordination 
identity 

pOO 

S{t) = / ip(t, t)T(t) dr, t> 0, (5.1) 

Jo 

with function ip(t, r) defined by 

1 /*7+ioo / \ 

<P(t,r) = — t > 0 (5 .2) 

2 7D Jj-ioo S 

where g{s) = g\(s) in case of problem 11.11 ) and g(s) = 52(a) in case of problem 
ED- The function ip(t,r) is a probability density function, i.e. it satisfies the 
properties 

/»00 

<p{t,T)> 0, / <p(t,T)dT = 1. (5.3) 

J 0 

Moreover, in the case of discrete distribution ED. there exists do G (0, 7t/2) such 
that (p{t,r) admits analytic extension to the sector |argi| < Qo and is bounded on 
each subsector \ argf| < 9 < 9q. 


Proof. According to (15.21) the Laplace transform of the function <p(t, r) with respect 
to t 





e st ip(t,r) dt, 


s, t > 0, 


is given by 

c p{s,r) = ^ e~ T9 ^ s \ s,r>0, (5.4) 

s 

Let T(t) be the solution operator of the classical Cauchy problem (11.71) and define 
an operator-valued function S(t) by (15.11) . Application of the Laplace transform 
gives by using (15.41) 


f oo noo r \ poo 

/ e~ st S(t) dt = / (p(s,t)T(t) dr = '■ - / e~ T9 ^T(r) dr = 

Jo Jo s Jo 


g{s) 


R(g(s),A). 


S 
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Comparing this result to (12.511 , it follows by the uniqueness of the Laplace transform 
that S(t) is exactly the solution operator of (11.111 if g = gi, resp. ( 11 . 21 ) if g = <? 2 - In 
this way ed is established. 

Let us prove now that <p(t,r) is a probability density function. Since g{s)/s £ 
CMJ- and e~ Tg ^ s ' > £ CMJ- as a composition of the completely monotone in x 
function e~ TX and the Bernstein function g(s), (15.41) implies that <f{s,r) £ CMT 
as a product of two completely monotone functions (see Proposition 2.1(a),(b) and 
Theorem 3.1(d),(e)). Then Bernstein’s theorem implies (p(t,r) > 0. 

The second identity in (15.31) can be proven in various ways. Here we show that 
it is a particular case of ED- Indeed, in the trivial case A = 0 the constant in t 
function u(t) = a satisfies equations ( 11 . 11 ) . ( 11 . 21 ) and (El) and by the uniqueness 
of the solution it follows that this is the solution in this case, i.e. S(t) =T(t ) = /. 
Then (15.11) reduces to the desired identity. 

In (16j, Section 4.2, it is proven that for g(s) = gi(s) in the case of continuous 
distribution the function (15.21) is well defined locally integrable function. In an 
analogous manner this can be done also for the continuous distribution variant of 
<72 (s). So, it remains to prove the last part of the theorem concerning the discrete 
distribution case. Define 

0o = nrin{(l/a — l) 7 r/ 2 , 7 r/ 2 } — e, (5-5) 

where £ > 0 is small enough, such that 0o > 0. According to m , Theorem 0.1, it 
suffices to prove that the function <J5(s,r) admits analytic extension to the sector 
I args| < n /2 + 9 0 and s^(s,t) is bounded on each subsector | args| < 7 r /2 + 6, 9 < 
9q. Theorem 3.1(f) states that g(s) can be extended analytically to £, r , thus this 
holds also for the function i)j(s,t). Take s such that | args| < tt/2 + 9 1 9 < 9q- 
Then, applying again Theorem 3.1(f), it follows 

| argg(s)| < a\ args| < 7t/2 — ae. 

Therefore, g(s) = for some p > 0, |</>| < 7 r /2 — ae, and thus 

|s£(s,r)| = \g{s)e~ T9(s) \ < pe~ Tpcos * < pe~ ap < (ea)" 1 , 
where a = rsinae > 0. With this the proof is completed. □ 

In the particular case of single term equation, such a subordination formula is 
given in [2], Theorem 3.1, where the function ip(t,r) is expressed in terms of a 
special function of Wright type. 

In [16j the subordination relation m is proven for the continuous distribution 
case of problem to with A being the Laplace operator. 

Note that the subordination identity ED shows again that in ordered Banach 
space positivity of T(t) implies positivity of S(t) . 

6 . Conclusion 

An approach is proposed for the study of the distributed order fractional evolu¬ 
tion equations in the Caputo and the Riemann-Liouville sense, rewriting them as 
abstract Volterra integral equations. The obtained results are based mainly on the 
properties of the kernels of these integral equations, and especially those related to 
complete monotonicity. 

The results proven in this paper hold for a large class of fractional evolution 
equations, involving several types of fractional differentiation of order less than one 
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as well as various possibilities for the operator A: e.g. the Laplace operator, general 
second order symmetric uniformly elliptic operators, operators leading to the so- 
called time-space fractional equations, such as: space-fractional derivatives (e.g. in 
the Riesz sense), fractional powers of the multi-dimensional Laplace operator, other 
forms of fractional Laplacian (see e.g. El]), fractional powers of more general 
elliptic operators, etc. 

The developed technique is also applicable to other related problems, for example 
the Rayleigh-Stokes problem for the generalized second grade fluid with fractional 
derivative model, see e.g. [3 m, or to more general abstract Volterra integral 
equations with kernel k(t), which Laplace transform k(s) is well-defined for s > 0 
and is such that ( k(s))~ 1 is a Bernstein function. 
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